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ABSTRACT 

Using the near-infrared (NIR) observations of the afterglows of the high redshift (5 < z < 25) gamma- 
ray bursts (GRBs) that wiU be detected by the Swift satehite, we discuss a way to study the cosmic 
reionization history. In principle the details of the cosmic reionization history are well imprinted in the 
NIR spectra of the GRB afterglows. However the spectroscopy with a space telescope is required to 
obtain such an information for a very high redshift (z > 15) unless the neutral fraction of the high- 
er universe is less than 10~^. The broad-band photometry has the higher sensitivity than that of the 
spectroscopy, so that the NIR photometric follow-up of the GRB afterglows is very promising to examine 
the cosmic reionization history. A few minutes exposure with a 8-m class ground-based telescope of the 
afterglow of the high-z GRBs will reveal how many times the reionization occurred in the universe. 

Subject headings: cosmology: observations — gamma rays: bursts — intergalactic medium — radiative 
transfer — techniques: photometric 



1. INTRODUCTION 

The evolutionary study of the intergalactic medium 
(IGM) is now one of the most active fields in astrophysics. 
Especially, the cosmic reionization history of the IGM 
attracts many researchers in recent years (see Loeb & 
Barkana 2001 for a review). Gunn-Peterson trough short- 
ward of the Lya resonance (Gunn & Peterson 1965) in the 
spectrum of quasars with redshift z ~ 6 indicates that the 
end of the reionization epoch is z '--^ 6 (Becker et al. 2001; 
White et al. 2003). On the other hand, the recent obser- 
vation of the polarization of the cosmic microwave back- 
ground (CMB) by WMAP suggests that the beginning of 
the reionization is z ^ 20 (Spergel et al. 2003; Kogut et 
al. 2003). One might think that there is a discrepancy 
between both observations. 

The detailed simulation of the six-dimensional radia- 
tive transfer shows that the reionization process proceeds 
slowly in an inhomogcneous universe (e.g., Nakamoto, 
Umemura, & Susa 2001). If the reionization began at 
z ~ 20, the neutral fraction in the universe decreased grad- 
ually, and the universe was ionized almost completely at 
z ~ 6, the apparent discrepancy can be resolved (Wyithe 
& Loeb 2003a; Haiman & Holder 2003; Ciardi, Ferrara, & 
White 2003; Chiu, Fan, & Ostricker 2003; Sokasian et al. 
2003; Onken & Miralda-Escude 2003). The scenario that 
the initial partially ionized epoch was simply followed by 
the complete ionization can explain the large Thomson 
scattering opacity observed from the CMB polarization 
and also is consistent with the previous picture in which 
the end of the reionization is z ~ 6. 

Cen (2003a,b) and Wyithe & Loeb (2003b) recently 
proposed a new scenario of the reionization; the uni- 
verse was reionized twice. Even in an inhomogcneous 
universe, enough strong intensity of the ultra-violet (UV) 
background radiation can ionize the universe quickly. At 
z ~ 20, the first reionization was made by the Population 



HI (Pop III) stars (metal-free stars) with a top-heavy ini- 
tial mass function (IMF) which yields a much higher UV 
emissivity than that of normal Pop I and II stars. Then, 
the universe was partially recombined at z ~ 15 when 
the transition from the Pop III to II was occurred and 
the UV emissivity was suddenly suppressed because of the 
different IMF. Finally, the UV photons from the Pop II 
stars increased gradually and ionized the universe again 
at z ~ 6. 

Gen's scenario can also resolve the discrepancy, but it is 
different from others in the history of the cosmic reioniza- 
tion; the first complete ionization was followed by the par- 
tially ionized epoch and then the second complete ioniza- 
tion. We should assess whether Gen's scenario is favored 
by using other observations. In this paper, an assessment 
by using the afterglows of the gamma-ray bursts (GRBs) 
is discussed. 

The usefulness of GRBs to investigate the high-z uni- 
verse is pointed out by many authors (e.g.. Lamb & Re- 
ichart 2000; Ciardi & Loeb 2000; Barkana & Loeb 2003). 
It is strongly suggested that the long duration GRBs arise 
from the collapse of a massive star (Galama et al. 1998; 
Ucmura et al. 2003; Hjorth et al. 2003; Price et al. 2003; 
Stanek et al. 2003; Matheson et al. 2003). Hence, the 
GRBs can occur at very high-z once the massive stars are 
formed. For example, the Swift satellite^ is expected to 
detect ~ 10 GRBs per year occurring at z > 10 (Lamb & 
Reichart 2000). Furthermore, if we fix an observing time 
from the prompt emission, the observed afterglow flux does 
not become so faint even for the extremely high-z GRBs 
because the earlier phase of the afterglow is observed in 
the cosmological rest frame (Ciardi & Loeb 2000). 

In this paper, we show that the near- infrared (NIR) pho- 
tometric follow-ups of the GRB afterglows are very useful 
to investigate how many times the cosmic reionization oc- 
curred. Although many techniques to prove the reioniza- 



^ Research Fellow of Japan Society for the Promotion of Science. 
^ http://swift.gsfc.nasa.gov/ 
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tion history have been proposed so far. for example, hydro- 
gen 21 cm hne tomography (Madau, Mciksin, & Rees 1997; 
Carilli, Gnedin, & Owen 2002; Furlanetto & Loeb 2002) ^ 
Lya damping wing measurement (Miralda-Escude 1998; 
Barkana & Loeb 2003), metal absorption lines (Oh 2002; 
Furlanetto & Loeb 2003), CMB polarization anisotropy 
(Haiman & Holder 2003), and dispersion measure in GRB 
radio afterglows (loka 2003; Inoue 2003), the NIR photo- 
metric colors may be the most promising technique using 
the current facilities. 

In order to investigate the reionization history using our 
photometric method, we need to determine redshifts of 
GRBs by other ways, e.g., detections of iron lines in X-ray 
afterglow spectra (Meszaros & Rees 2003) and of the Lya 
break in NIR spectra, or empirical methods by using only 
the 7-ray data (Fenimore & Ramirez-Ruiz 2000; Norris et 
al. 2000; loka & Nakamura 2001; Amati et al. 2002; At- 
teia 2003; Murakami et al. 2003; Yonetoku et al. 2003). 
Even if redshifts of GRBs arc unknown prior to follow-up 
observations, it is worth performing NIR photometry as 
early phase as possible since ^ 10% of GRBs are expected 
to be located at 2; > 10 (Bromm & Loeb 2002). In practice, 
such follow-up observations against every GRB are possi- 
ble. Even after the early NIR follow-up, we will be able to 
determine the redshifts. In this paper, we show that the 
near-infrared (NIR) photometric follow-ups of the GRB af- 
terglows provide a significantly important information on 
the cosmic reionization history. 

The structure of this paper is as follows: we start from 
a brief summary of the IGM opacity by the neutral hy- 
drogen in section 2. In section 3, the possible scenario 
of a very low neutral fraction at very high-z is discussed. 
Then, we examine the NIR spectra and colors of the GRB 
afterglows in section 4 and 5, respectively. Finally, we dis- 
cuss a way to obtain the cosmic reionization history and 
advantage/disadvantage of our method in section 6. 

We adopt a standard set of the ACDM cosmology 
throughout the paper: Hq = 70 km Mpc"-^, CIm = 0.3, 
^Ia = 0.7, and Q), = 0.04. 

2. IGM OPACITY 

Suppose an observer at z = who observes a source 
with z = zs at the observer's frequency vq. The radiation 
from the source is absorbed by the Lyman series lines and 
the photoionization of the intervening neutral hydrogen 
(Gunn & Peterson 1965). The hydrogen cross section at 
the rest frame frequency u is 



and 



£7Hl(!^) = (TLc(l^) + , 
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where ctlc('^) is the cross section for the Lyman contin- 
uum photons {hi' > 13.6 eV) and <Ti{i') is the cross section 
for the i-th line of the Lyman series, i.e., i = Lya, Ly/3, 
Ly7, etc. Here we consider the lines up to i = 40. The 
optical depth for the observer's frequency vq is given by 
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^ In Figure 9 of Cen (2003a), the neutral fraction decreases up to 10~*. 
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where nm{z) is the number density of the neutral hydro- 
gen at the redshift z, H{z) is the Hubble constant at the 
redshift z, and c is the light speed. 

The cross section for the Lyman continuum photons is 
given by <TLc(t^) = o'l('^/^'l)~^, where ctl = 6.30 x 10"-^* 
cm^ and the Lyman limit frequency = 3.29 x 10^^ Hz 
(Osterbrock 1989). Thus, we obtain 
3 
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where Xhi{z) is the neutral fraction at the redshift z and 
Zmin = niax[0, (i'l/i'o) — 1]- The neutral fraction xm{z) 
is defined as nHi(z)/nH(z); where nH(z) = nH,o(l + z)^ 
is the cosmic mean number density of hydrogen atom at 
the redshift z with nn.o being the present one. The effect 
of the inhomogeneity of the universe can be expressed in 
the definition of xhi (see section 3 below). Since we need 
the condition Zmin < -^s to integrate Eq. (5), tI'^{zs) = 
when uo < y-L/{l + 2s)- 

The line cross section (Ji{v) has a very sharp peak at 
the line center frequency. The width is well character- 
ized by the Doppler width, Ai^D- That is, ai{v) ~ ai for 
\v — < Az/D, and ai{v) ~ for \v — Vi\ > Ai/d, where 
CTj and Ui are the cross section and the frequency of the 
t-th line center, respectively. Since v = uq{\-\- z),we inte- 
grate the rhs of equation (4) over a narrow redshift width 



of A(l + z) = 2(Ai/D/t'i)(l + Zi) around 1 
Therefore, we obtain 
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where we have substituted 
which e is the electric charge, fi is the absorption oscillator 
strength of the i-th line, and nie is the electron mass. The 
values of fi and Vi are taken from Wiese, Smith, & Glennon 
(1966). This time also we are restricted by Zi < zs- Thus, 
the above equation is valid only when vq > i/j/(l -|- zs), 
otherwise t^^(zs) = 0. 

For the Lya line, the opacity becomes 

3/2 

(7) 
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where we approximate H{z) « HoQlp{l + z)^/^ and 
adopt nHi(z) = .XHi(2)nH,o(l + z)^. This is the Gunn- 
Pcterson optical depth (Gunn & Peterson 1965; Peebles 
1993). From this equation, we realize that if xhi ^ 10~^, 
the IGM is opaque for the photon bluer than the source 
Lya line. Conversely, we can estimate Xui if t^^" is de- 
termined observationally (e.g., Becker et al. 2001; White 
et al. 2003). 

3. NEUTRAL FRACTION AND STAR FORMATION RATE 

In the following sections, we will assume a neutral frac- 
tion of hydrogen, xm ~ 10~^ at z ~ 20 for calculations 
of spectra and colors of GRB afterglows. Before such cal- 
culations, we show, here, that xui ~ 10~^ is possible at a 
very high-z universe.^ 

However, it is an artificial limit in his calculation (Cen 2003, private 
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In general, the neutral fraction, xmiz), is spatially vari- 
able. However, we discuss only its mean value for simplic- 
ity. The averaged Xui is determined by the UV intensity 
of the background radiation and the inhomogcncity of the 
hydrogen number density at the rcdshift z. Even at very 
high-z universe, the recombination time-scale is much less 
than the Hubble time-scale. Thus, we assume the ioniza- 
tion equilibrium. In that case, the mean neutral fraction 
for a highly ionized medium (i.e. Xm ^ 1) is given by 
.thi ~ Cnna/Tm, where C = (n^)/(riH)^ is the clumping 
factor, riH is the number density of hydrogen nuclei, a is 
the recombination coefficient, and Fhi = / ai,c{i^)cnvdi' is 
the HI photoionization rate where is the photon num- 
ber density per frequency. Therefore, the UV background 
intensity is required to estimate xm even if we know the 
clumping factor from the model of the cosmological struc- 
ture formation. 

How much photons are required to control the hydrogen 
neutral fraction in a very low level? Approximately, the 
HI photoionization rate is estimated to be Fhi ^ o"Lcnion 
with nion being the number density of ionizing photons. 
Hence, we find 

C('^\. (8) 



riH (Jhcxm \ xhi 
where the case B recombination rate a = 2.73 x 10"^"^ cm^ 
s"-'' is adopted (Osterbrock 1989). Now. we are interested 
in a very high redshift universe {z ~ 20). In there, the 
clumping factor is an order of unity. Therefore, we find 
that only one photon per a hydrogen nuclei is suflacient to 
keep .XHi 10-*^. 

Then, we examine how much stars are required to main- 
tain the photon density. Since a proton recombines with 
an electron in a certain time-scale, continuous supply of 
ionizing photons is needed to keep rijon ^ nn. Since the 
recombination time-scale is ~ l/nna for xm <C 1, the 
required photon emissivity per unit volume is estimated 
to be Cion ^ '^ion'^^H'^ ^ Cn^Q.'^ / xniCTi^c. On the other 
hand, the emissivity is given by eion = ./osceLC/fSFR(l + z)^, 
where /esc is the escape fraction of Lyman continuum from 
primordial galaxies, clc is the Lyman continuum photon 
emissivity per unit stellar mass, and psFR is the star for- 
mation density per unit time per unit comoving volume. 
Therefore, the required star formation density is 

PSFR 



Mq yr-i Mpc-3 ^'^'^ \f, 



10- 



0.1 \ /lO'^iphM©-^ 



1 



, - (9) 

Although the photon emissivity at high-z is very uncer- 
tain because we don't know the stellar mass distribution, 
we estimate the emissivity by Starburst 99 model (Lei- 
thercr et al. 1999)'' by assuming the Salpeter mass func- 
tion with various mass ranges. The estimated values of 
clc are summarized in Table 1. The Pop HI stars are 
likely to have a top-heavy mass function (e.g., Nakamura 
& Umemura 2001). Hence, the case of 10 100 Mq in Table 
1 may be suitable. In this case, the required star forma- 
tion density to maintain xui ~ 10~^ is ~ 0.05 Mq yr~^ 
Mpc~^ (comoving) when f^sc = 0.1 and C = 1. 

Is this star formation density possible? The latest ob- 
servations suggest that the star formation density retains 

^ http://www.stsci.edu/science/starburst99/ 



a level of ^ 0.1 Mq yr~^ Mpc""^ from z ~ 1 to z ~ 6 
(Giavalisco ct al. 2003). Such a level of star formation 
may be kept toward more high-z universe. Moreover, a 
semianalytic model shows 0.01-0.1 Mq yr~^ Mpc~^ at 
z ^ 20 depending on the assumed star formation efficiency 
(Somerville & Livio 2003). Thus, we can sufficiently ex- 
pect ~ 0.05 Mq yr-^ Mpc~^, and then, xm ~ 10~^ at 
z ~ 20. In addition, we note that the escape fraction may 
be much larger than 0.1 assumed above if the Pop HI stars 
are formed in low-mass halos. 



4. NEAR INFRARED SPECTRA 

Let us discuss the observed spectra of the GRB after- 
glows in the NIR bands. To do so, an afterglow spectral 
model is required. We adopt a simple afterglow model; the 
synchrotron radiation from the relativistic shock (Sari, Pi- 
ran, & Narayan 1998). More specifically, we adopt equa- 
tions (l)-(5) in Ciardi & Loeb (2000) who take into ac- 
count the effect of the cosmological redshift. The adopted 
parameters are the magnetic energy fraction of cb = 0.1, 
the electron energy fraction of = 0.2, the spherical shock 
energy oi E = 10^^ erg, the ambient gas number density 
of n = 10 cm^"^, and the power-law index of the electron 
energy distribution p = 2.5. 

First we consider the observed afterglow spectra in the 
hypothetical perfectly neutral universe for comparison. In 
Figure 1, we show the expected afterglow spectra in the 
neutral universe observed 1 hour after the burst in the 
observer's frame. Due to the Lya line absorption, the 
continuum bluer than the Lya line in the source frame 
(the observed wavelength 0.1216[1 + zs] M™) is completely 
damped. Thus, we can find the Lya break clearly. From 
the observed wavelength of the Lya break, we can deter- 
mine the redshift of the GRBs. If we observe the after- 
glow through a filter, the radiation from the source with 
the redshift higher than the characteristic redshift of the 
filter cannot be detected because of the Lya break. For 
example, the efi'ect starts from — 8 for the J-band and 
the source with zs > 11 cannot be seen through the filter, 
i.e. J drop-out. These characteristic redshifts arc sum- 
marized in Table 2. However, we have a chance to see 
the source beyond the drop-out redshift if the universe is 
highly ionized as we will show later. 

Next we examine what is observed if the very high-z uni- 
verse is ionized completely as proposed by Ccn (2003a, b). 
Let us set the neutral fraction, xhi, of the universe in the 
redshift range 15 < z < 20 to be very small and xhi ^ 1 
for z < 15 and z > 20. That is, we assume that the Pop 
HI stars ionized the universe at z = 20 and the sudden 
change of the IMF due to the transition from the Pop HI 
to II was occurred at z = 15. The neutral fraction is de- 
termined by the background UV intensity produced by the 
Pop HI stars. However, the intensity is quite uncertain, 
so that we choose two cases as xhi = 10^^' and 10^^ for 
instance (see also Fig. 5 in section 5). We will present a 
way to constrain xm from the observations later. 

In Figure 2, we show the expected spectra of the GRB 
afterglows. The solid, short-dashed, long-dashed, dot- 
dashed, and dotted curves are the expected afterglow spec- 
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tra of the GRB at zg = 22, 20, 18, 15, and 13, respec- 
tively. The observing time is assumed to be 1 hour af- 
ter the prompt emission in the observer's frame. We find 
the clear Lya breaks in the spectra. Since xhi = 10~^ 
(or 10~^) at 15 < z < 20 is assumed in the panel (a) 
(or (b)), the IGM opacity is an order of unity (or 0.1) in 
the redshift range (see eq.[7]). Thus, the continuum bluer 
than the Lya break of the GRBs with zs > 15 still re- 
mains of the order of 10 /xJy (or 100 /xJy) in ~ 2-2.5 /xm, 
i.e., in K-hand (sec the thin solid curve indicated as K; 
we adopt the filter system of Bessell & Brett 1988). The 
spectral break at 1.94(=0.1216[1-|-15]) /xm is due to the 
neutral hydrogen below z = 15. Thus, the continuum less 
than the break wavelength in the observer's frame from the 
Zs > 15 source is completely extinguished. The spectra of 
the source with < 15 are the same as those shown in 
Figure 1. For zs = 22 case, the spectrum shows the second 
break at 2.36(=0.1026[l+22]) ^m. This is the Ly/3 break 
due to the neutral hydrogen near the GRB. Here we have 
assumed .thi ^ 1 for z > 20. 

The structure corresponding to the rcionization history 
appears in the spectra of the GRB afterglows as shown 
in Figure 2 (sec also Haiman & Locb 1999). If we could 
detect the continuum rising at 2.55(=0.1216[l-|-20]) /um in 
the spectrum of the GRB with zs = 22, we would find 
the starting epoch of the first rcionization as z = 20. On 
the other hand, the end of the first rcionization, in other 
words, the transition epoch from the Pop III to II is real- 
ized from the spectral break at 1.94(=0.1216[H-15]) /um. 

The suitable band to determine the ionization history 
depends on the redshift of the rcionization epoch. From 
Figure 1 (see also Tabic 2), wc find that the /, J, H, 
and K-hand spectroscopies are suitable for the rcioniza- 
tion epoch at z ~ 5-7, 8-11, 11-14, and 15-20, respec- 
tively. In any case, observations to detect the spectral 
signatures in the NIR afterglow spectra of the GRBs are 
strongly encouraged. If Xm is smaller than 10~^, we can 
clearly see the difference from Fig. 1 and confirm the dou- 
ble rcionization. 

In the above discussion, we assumed xhi = 10~^ and 
10~^. Let us argue what will happen for different val- 
ues of Xhi- If a;Hi ^ 10~^, the remaining flux decreases 
exponentially because the IGM opacity becomes much 
larger than unity; for example, the flux is about 1 nJy 
for a;Hi = 6 x 10~^, and about 1 pJy for xhi = 10~^. It is 
quite diflicult to detect such a low-level flux. 

5. NEAR INFRARED COLORS 

The photometric observations are available easier than 
the spectroscopy. We examine the expected apparent NIR 
colors of the GRB afterglows. Although we can discuss the 
apparent magnitude of the afterglows in one photometric 
filter, their dispersion is very large because the luminosity 
of the afterglows depends on many uncertain parameters 
siich as the jet opening angle, the ambient matter den- 
sity, the magnetic energy fraction and the relativistic elec- 
tron energy fraction. On the other hand, the dispersion of 
the apparent colors can be quite small because the color 
does not depend on the absolute luminosity but only the 
spectral shape which does not change significantly in the 
observed NIR bands. 

^ All magnitudes in this paper are the Vega system. 



The apparent magnitude^ in a filter band denoted as i 
is defined by 

mi = -2.51ogFi/Fi,o, (10) 

where Fi^ is the zero point flux density of the filter and 
Fi is the mean flux density through the filter which is 



Fi = 



(11) 



where Tj^^ is the transmission efficiency of the filter and 
/ye"^" is the fiux density entering the filter. If there is no 
intervening absorber between the source and the telescope, 

= 0. 

The observed color between two filter bands, i and j 
(the filter i is bluer than the filter j), is given by 



rrii — nij = (mj — mj)™^ + {Ami — Amj) , 



(12) 



where (mj — mj)™^ is the intrinsic color of the source, and 
Anii and Amj are the absorption amounts in the filters 
i and j, respectively. When we consider the NIR filter 
bands and high-z GRBs (for example z = 15), the intrin- 
sic afterglow spectrum is predicted to be proportional to 
J/-1/2 from ~ 1 minute to several hours after the burst 
occurrence and proportional to i/"''/^ for later time in the 
standard afterglow model (Sari, Piran, & Narayan 1998; 
Ciardi & Loeb 2000). Other parameters adopted are de- 
scribed in the beginning of the section 3. In Table 3, we 
tabulate the intrinsic colors of the sources for the two cases 
of the spectral shape. 

The absorption amount in the filter i is 



Ami = mi- mf* = -2.5 los 



/ Ti^ufi, 



'dv 



J T.^ufudv 



(13) 



where mJ"* = mi{T„ = 0) is the intrinsic (no absorp- 
tion) apparent magnitude. In a certain band width Ai>, 
the difference in the optical depth At is estimated as 
\At/t\ = 3/2(1 + z)-^\A{l + z)| - (1 z)-^\Aiy/iy\ from 
equation (7). Since Az/ of the filter transmission is smaller 
than the effective frequency of the filter, i.e. Avjv < 1, 
and also z > 5, At/t <C 1, that is, the term e^'^'^ in 
the integral of the numerator in equation (13) can be re- 
garded as almost constant. Hence, we obtain approxi- 
mately Ami ~ 1.086Teff, where Tefj is the effective IGM 
opacity in the filter i. 

Now we can estimate the observed color by equation (12) 
if Ami is known. To know Ami is equivalent to know the 
IGM effective opacity Tes- This opacity is one between the 
redshift at which the Lya break comes into the filter band 
width (zLy^ i„, see Table 2) and the source redshift (zs) be- 
cause the neutral hydrogen in z < -2Lya,in cannot absorb 
the photons passing through the filter i. We note that the 
neutral hydrogen lying beyond the redshift at which the 
Lya break goes out of the filter band width (z^y^, see 
Table 2) absorb the photons through the filter i because 
of the higher-order Lyman series lines like Ly/3, Ly7, etc., 
and the photoionization process. As a result, is de- 
termined by the neutral fraction a;Hi in the redshift range 
zhya,in ^ ^ < zs- Sincc o^Hi at high-z is uncertain, we as- 
sume that Xhi is constant in the above range for simplicity. 
The real xm might vary significantly in the redshift range, 
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so that the assumed xm should be regarded as an effective 
mean value including such a variation (hereafter x^). 

In Figures 3a-3d, we show Am, for the /, J, H, and 
-fT-bands as a function of the source redshift. The contin- 
uum in the observer's i-band is not absorbed by the IGM 
neutral hydrogen at all when the source redshift is less 
than about 25. Although we assumed that the spectral 
shape is proportional to u~^/^ in the panels, the results 
are much robust for the change of the spectral shape as 
noted above. The solid curves in these panels are loci of 
Ami for a given as a function of the source redshift 
zs- For example, since ^Lya.in = 5 for 7-band, the IGM ab- 
sorption in /-band is about 5 mag for the source at zs = 7 
if the effective neutral fraction in the redshift range 
5 < < 7 is 10-^ 

Two dotted vertical lines in each panel of Figure 3 in- 
dicate the redshifts when the Lya break enters and goes 
out of the each band width (zLy^^in and ZLyQ,_outJ respec- 
tively). As seen in Figure 1, if the IGM is significantly 
neutral, the source with zs > ^Lya^out cannot be detected 
through the filter i (i.e. drop-out). However, we can de- 
tect such a source if the universe is highly ionized. This 
is because the continuum below the Lya break remains as 
shown in Figure 2. 

Suppose we observe a source with zs > ^Lya out through 
the i and j filters and assume that the radiation through 
j filter is not affected by any intervening absorption (i.e. 
Arrij = 0). Prom equation (12), the expected magnitude 
through the i filter is 

rrii = rrij + {rrii — mj)™' -I- Ami ■ (14) 

For example, we consider the case of i = /, j = L, and 
Zs = 7. We find that the apparent L magnitude of the af- 
terglow of the GRB at zs = 7 for 1 hour after the prompt 
burst is about 14 mag from Figure 4. Thus, the apparent 
/ magnitude is expected to be 22 mag because the intrin- 
sic / — L = 3.1 mag for 1 hour (i.e. the case oc in 
Table 3) and A/ ~ 5 mag if x'^^ = 10"^ in the redshift 
range 5 < z < 7. We can reach 5-(t detection of the source 
with I = 22 mag by only three minutes exposure with a 
8-m class telescope. Interestingly the assumed is sim- 
ilar to the value reported by White et al. (2003) from the 
Gunn-Peterson trough in the spectra of z > 6 quasars. 

Similar argument can be done for other bands. There- 
fore, in general, the detection of the source with zs > 
-^Lya.out through the filter i is the very good evidence that 
the universe in the redshift range -^Lya.in ~ ~ -^^s is highly 
ionized, i.e. <C 1 in that redshift range. Conversely, 
the detection enables us to estimate Tes and a;|^. Finally, 
we note here that Figure 3 is also useful for any other 
sources (e.g., QSOs) because Ami is almost independent 
of the source spectrum. 

6. DISCUSSIONS 

6.1. Was the universe reionized twice? 

Here we discuss how we can confirm or refute Gen's sce- 
nario: the universe was reionized twice. In the scenario, 
the first complete ionization at z ^ 20 is followed by the 
partially ionized epoch at z ~ 10. Therefore, we should 
check whether the neutral fraction at z ~ 20 is very low 
or not and whether the fraction at z 10 is almost unity 



or not. To do so, the best observation is the spectroscopy 
in the NIR bands. As shown in Fig. 2 of section 3, the 
reionization history is imprinted in the observed contin- 
uum. However, the sensitivity of the spectroscopy is in 
general much less than that of the photometric observa- 
tions. Hence we discuss the way using the NIR photome- 
tries. 

In the previous section, we have shown that the detec- 
tion of the GRB afterglows through a filter i beyond the 
Lya drop-out redshift (zl^y^^^^^.) proves the ionization of 
the universe around z^^^ ^j^^ . From the characteristic red- 
shifts for the NIR filters summarized in Table 2. the /, 
J, H, and K filters are the suitable to check the ioniza- 
tion state at z ~ 5-8, 8-11, 11-15, and 15-20, respec- 
tively. Thus, the null detection of the GRB afterglows of 
zg ^ 11 in J-band indicates an high neutral fraction in 
8 < z < 11. On the other hand, we detect the GRB after- 
glows of zs ~ 20 in K-h&nd if the IGM in 15 < z < 20 is 
ionized. / and i?-band surveys are also very important to 
assess the reionization history of the universe. We can con- 
strain the latest reionization epoch by observing the GRB 
afterglows at z > 6 through /-band. In summary, we can 
examine the reionization history by checking whether the 
high-z GRB afterglows drop out of the NIR filters or not. 

In the rest of this subsection, we discuss what is the 
difference between Ccn's scenario and others. To demon- 
strate the main feature, we assume two schematic reion- 
ization histories; (1) single gradual reionization and (2) 
double sudden reionizations, which are shown in Figure 5 
as the solid and dashed curves, respectively. These his- 
tories are based on two observational constraints; (1) the 
neutral fraction of hydrogen xui ~ 10~^ at z ~ 6 from the 
Gunn-Peterson trough found in the z > 6 quasars spectra 
(Becker ct al. 2001; White et al. 2003), and (2) the begin- 
ning of the reionization is z ~ 20 from the large opacity 
of the electron scattering suggested by WMAP (Kogut et 
al. 2003). For the double reionizations, we also consider 
different values of xui in the first reionization epoch. 

In Figure 6, we show the expected NIR colors as a 
function of the source redshift for the afterglow spectrum 
oc v~^/^ case (the observing time less than several 
hours, see Table 3). We find differences between the sin- 
gle and double reionizations in the I — J (panel [a]) and 
K ~L (panel [d]) colors. In panel (a), the GRB afterglows 
up to zs ~ 8 can be seen in both of the / and J bands 
for the single reionization case, whereas the I — J color 
of Zg > 7.1 afterglows diverges in the double reionization 
case, i.e. the sources drop out of the /-band. This re- 
flects the difference of the increasing rate of the neutral 
fraction around z ~ 6 in two reionization histories. The 
drop-out redshift in the double reionization case is deter- 
mined by A/3(H-zs,drop) = AQ,(1 + Zi.oion), where Aq and A^^ 
are the rest-frame wavelength of the Lya and Ly/j lines, 
and Zrcion is the sudden reionization redshift (Haiman & 
Loeb 1999). In our case, Zreion = 6. It is worth to noting 
that the GRB afterglows beyond the drop-out redshift of 
the /-band (zLy„ out = 6.6) can be seen through the filter 
in both of reionization histories. 

In panel (d) of Figure 6, we find a significant difference 
between the two reionization scenarios. While the after- 
glows beyond the drop-out redshift of the A'-band really 
drop out of the filter for the single reionization, we can 
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see such afterglows in the twice reionized universe. On the 
other hand, the afterglows with zs > 23 drop out of the 
i^-band even in the double reionization case because the 
Ly/3 break goes out of the filter transmission width, i.e. 
A/3(l + 2s,drop) = A^a^h A„(l + 2;/^y„_out)]' where A^^^ is 
the maximum wavelength of the if-band filter. We note 
here that 2;reion(= 20) > ^^Lya.out ™ time, whereas 
•^■reion (= 6)< 

•^Lya,out above case. 

Any difference between the single and double reioniza- 
tions cannot be found in the J — H (panel [b]) and the 
H — K (panel [c]) colors because the neutral fractions in 
both cases are high enough (> 10^"') to extinguish the 
continuum bluer than the Lya break completely. That 
is, the GRB afterglows beyond the drop-out redshift can- 
not be seen in the J and i?-bands for both of reionization 
histories. 

Let us summarize how to confirm or refute Cen's sce- 
nario: We can conclude that the universe was reion- 
ized twice if (1) the GRB afterglows with zs > z^ya.ont 
is detected in the iiT-band and (2) the afterglows with 
> ^Lya.out or Zs > z^ya^ont ^rop out of the J or H- 
bands. However, the null detection of the zg > ZLya out 
GRB afterglows in the X-band does not reject Cen's sce- 
nario at once. The null detection only shows the neu- 
tral fraction at z ~ 20 larger than ~ 10^''. In any case, 
the deep and prompt /T-band photometry of the high-z 
GRB afterglows is useful to examine the ionization state 
at z ~ 20 very much. 

6.2. Comment on Lyman break technique 

In the above discussions, we assumed that the GRB red- 
shifts are known from other methods. Although the spec- 
troscopy of the optical afterglow or the host galaxies is 
used to determine the redshift of the low-z GRBs, it may 
be difficult for the high-z GRBs. The Fe line in the X-ray 
afterglow (Mcszaros & Roes 2003) or some empirical ways 
using the 7-ray data alone (e.g., Yonetoku et al. 2003) can 
be useful for the high-z GRBs. In addition, the search for 
the Lya break is considered to be useful. Indeed, the spec- 
troscopic detection of the sharp Lya cut-off (see Figs. 2 and 
3) is an accurate method to determine the redshift. How- 
ever, the color selection technique which is often used to 
find the z > 3 galaxies (Madau 1995; Steidel et al. 1996) 
may not be good for the high-z GRBs. This is because the 
GRB afterglows beyond the drop-out redshift of the con- 
sidered filter can be detectable if the universe is ionized 
enough. 

For example, if we use / — J > 5 mag as a criterion to 
select the high-z GRBs, the selected objects have z > 6.5 
indeed (see Fig. 3a and Table 3). However, a number of real 
z > 6.5 objects escape the criterion if the neutral fraction 
is less than ~ 10~^. Therefore, we cannot use a simple 
color selection technique for the high-z GRBs. Spectro- 
scopic observations to detect the Lya break feature , Fc 
lino (Mcszaros & Recs 2003) or empirical methods by us- 
ing only the 7-ray data (Fcnimorc & Ramirez-Ruiz 2000: 
Norris et al. 2000; loka & Nakamura 2001; Amati et al. 
2002; Atteia 2003; Murakami et al. 2003; Yonetoku et al. 
2003) are required. 

^ http;//subarutelescope.org/ 

^ http://astro.estec.esa.nl/SA-general/Projects/Planck/ 



6.3. Advantage and disadvantage of NIR colors method 

The largest advantage of the NIR color method is its eas- 
ier availability and higher sensitivity than other methods. 
The limiting magnitude with a high signal-to-noise ratio 
of the broad-band NIR photometry reaches much deeper 
than 20 mag by only a few minute exposure with a 8-m 
class ground based telescope. For example, the IRCS (In- 
frared Camera and Spectrograph) equipped with the Sub- 
aru telescope (Japanese 8-m class telescope^) can reach the 
2-fT upper limit of 22 mag by only five minutes exposure 
in X-band. This magnitude corresponds to AK ~ 5 mag 
based on the expected L magnitude of ~ 16 mag for the 
z ~ 20 GRB afterglows at 1 hour after the burst (Fig. 4) 
and K — L = 1.1 (Table 3), and also corresponds to the 
neutral fraction of 3 x 10~°. Thus, we can put this value 
as a lower limit of the neutral fraction at z ^ 20 for null 
detection. 

A more strict lower limit can be obtained if we detect 
the emission from the reverse shock (Sari & Piran 1999a, b; 
Gou et al. 2003) . Some GRBs show a very bright early af- 
terglow from the reverse shock. Interestingly, the apparent 
magnitude becomes ~ 5-6 mag brighter than that shown 
in Figure 4. In this case, we can reach AK ~ 10 corre- 
sponding to Xhi — 6 X 10~^. A very early observation is 
highly desired. 

From the detected magnitude or the limiting magnitude 
for the null detection, the neutral fraction or its lower 
limit in the corresponding redshift range can be estimated. 
However, the uncertainty of the obtained neutral fraction 
may be large if we have only one photometric data. This 
is because the apparent dispersion of the afterglow lumi- 
nosities is large. Even in that case, the uncertainty can be 
controlled in a low level if we use more than two photomet- 
ric data, i.e. colors, which arc independent of the absolute 
luminosity. This point is one of the important advantages 
of the NIR color method. Therefore, the NIR multi-colors 
follow-ups of the GRB afterglows are strongly encouraged. 

Miralda-Escude (1998) and Barkana & Loeb (2003) 
show that the detailed spectroscopy of the red damp- 
ing wing of the Lya break provides us with the optical 
depth for the line, i.e., the neutral hydrogen column den- 
sity to the source. However, the spectral resolution re- 
quired is A/AA ~ 5000(10^/rLya)- For a lower opacity, a 
much higher resolving power is needed. The limiting mag- 
nitude for such observations becomes significantly shal- 
low. That is, the method of damping wing measurement 
do not have sensitivity for a low opacity case. On the 
other hand, the NIR colors method is practically sensitive 
against T-^ya ^ 10- Therefore, these two methods are com- 
plementary each other; if xhi ~ 10~^ ('^Lya ~ 1), the NIR 
colors method becomes very useful, whereas the damping 
wing method is promising when .thi ^ 0.1 [r-Lya ^ 10^). 

Measurement of the CMB polarization anisotropy is also 
useful. Haiman & Holder (2003) show that we can distin- 
guish how much times the universe is reionized by using 
EE spectrum with Planck sensitivity in 3cr level. Since 
the sensitivity of WMAP is not enough, we must await 
the launch of Planck^ to reveal the reionization history by 
the CMB measurement. 
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Moreover, metal absorption lines like O i A1302 can be 
useful to prove the reionization history (Oh 2002; Furlan- 
etto & Loeb 2003). However, the expected equivalent 
width is very small as < 5 A. We require the spectroscopy 
with a resolving power ^ 5000. Unfortunately, observa- 
tions with a ground-based telescope may be difficult be- 
cause a huge number of Earth's atmospheric OH lines con- 
ceal the metal lines. Only JWST (James Webb Space Tele- 
scope)^ will have such a high spectral resolution among 
the future space telescopes having the NIR spectrograph. 
Thus, we need to await its launch to discuss the metal 
absorption lines. 

Technique using the hydrogen 21 cm absorption line is 
also proposed (e.g., Madau et al. 1997). However, the 
brightness of GRB afterglows is too faint to use them 
as a background light source for the present and future 
radio facilities because the absorption line is very weak 
(Furlanetto & Loeb 2002). Thus, we need to look for other 
sources. Although Carilli et al. (2002) proposed luminous 
high-^ radio-loud quasars as a candidate of the background 
source, it is very uncertain that they exist at z ^ 20. 

The dispersion measure in GRB radio afterglows may be 
promising in near future (loka 2003; Inoue 2003). If we ob- 
serve the radio afterglow at about 100 MHz within about 
1000 s after the burst occurrence, the delay of the arrival 
time of the low frequency photons may be detectable by 
the Square Kilometer Array. As well as the measurement 
of the red damping wing of the Lya line, this technique is 
sensitive to Xui ^ 0.1 (i.e., TLya ^ 10^). Thus, this tech- 
nique and our NIR color method are also complementary 
each other. 

The disadvantage of the NIR color method is the coarse 



redshift resolution. As discussed in section 4, the neutral 
fraction obtained by the color method is an averaged one 
in the redshift range of -JLya.in < z < zs- Thus, we cannot 
determine the so-called reionization redshift. Only the red- 
shift range in which the reionization occurred is obtained. 
However, even such a rough redshift resolution is enough 
to discuss whether the universe was reionized once, twice, 
or more. 

To know the detailed history of the reionization, the 
spe(;troscopy is needed. For this purpose, the Japanese as- 
tronomical satellite ASTRO-F^ can be useful. The satel- 
lite will have the spectroscopic sensitivity of ^ 30 /uJy 
around ^ 2 /im, which can detect the GRB afterglows of 
zs ~ 20 if xm - 10-^ (sec Fig.2a). SIRTF (Spccc In- 
fraRed Telescope Facility) may not be useful because it 
has the sensitivity only in the wavelength longer than 5 
/Ltm. In future, JWST is much promising. Finally, we note 
that the narrow-band photometry may be useful because 
its narrow transmission width provides us with a moder- 
ate resolution of redshift keeping a higher sensitivity than 
that of spectroscopy. 
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Fig. 1. Afterglow spectra in the hypothetical neutral universe. The solid curves are the afterglow spectra of the gamma-ray bursts at the 
redshift zg = 4-23 as indicated in the panel. The observing time is set to be 1 hour after the burst for all curves. The dotted curves are the 
filter transmissions of /, J, H, K, and L bands (Bessell & Brett 1988; Bessell 1990). These transmission curves are scale free. 
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Fk;. 2. — Expected near-infrared afterglow spectra of very high redshift gamma-ray bursts. The solid, short-dashed, long-dashed, dot- 
dashed, and dotted curves are the cases of the source redshift = 22, 20, 18, 15, and 13 respectively. The observing time is set to be 1 hour 
after the burst for all curves. The neutral fraction in 15 < z < 20 is assumed to be 10~® and 10^^ for the panels (a) and (b), respectively. 
The thin solid curve is the filter transmission of X-band (Bessell & Brett 1988). This transmission curve is scale free. 
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Fig. 3. — The intcrgalactic absorption in ncar-infrarcd bands. The vertical axis in cacli panel means the absorption amount in each filter 
(sec cq. [11]); (a) 7-band, (b) J- band, (c) i/-band, and (d) X-band. The solid curves in the panels are loci of the absorption amount as a 
function of the redshift of the source zs- The neutral hydrogen fraction, a^jjii is assumed to be constant between the redshift when the Lya 
break enters the filter transmission and zg- The assumed xhi is indicated on each curve. Two dotted vertical straight lines in each panel 
indicate the source redshifts at which the Lya break enters and goes out of the filter. 
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Fig. 4. — The expected apparent magnitude in L-band of the afterglows of the gamma-ray bursts (GRBs) as a function of the source 
redshift. The dashed, solid, and dotted curves are the cases of 1 minute, 1 hour, and 1 day after the burst occurrence in the observer's frame, 
respectively. The spectral model of the GRB afterglows by Ciardi & Loeb (2000) is adopted. The assumed parameter set is described in the 
first paragraph of section 3. 
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Fig. 5. — Examples of the reionization history. The sohd curve is a single gradual reionization case. The dashed curve is a double 
reionizations case as suggested by Cen (2003a). 
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Fig. 6. Expected near-infrared colors of gamma-ray burst afterglows, (a) / — ,7 colors, (b) J — H colors, (c) H — K colors, and (d) K — L 
colors. The solid and dashed curves correspond to the cases of the single and double reionizations depicted in Figure 5, respectively. The 
solid and dashed curves in panels (b) and (c) are perfectly superposed. In panel (d), some cases of the neutral hydrogen fraction in the first 
reionization are shown. Two dotted vertical straight lines in each panel indicate the source redshifts at which the Lya breaJj enters and goes 
out of the indicated filter. 
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Table 1 

Star formation density required to maintain xhi ~ 10 



Stellar mass range 




PSFR 


(Mq) 


(lO^'i ph Mq-I) 


(M0 yr-1 Mpc-3) 


10-100 


2 


0.05 


1-100 


0.7 


0.1 


0.1-100 


0.3 


0.3 



Note. — We calculate the Lyman continuum photon emissivity 
per unit stellar mass, clc, by Starburst 99 model (Lcitherer et 
al. 1999) by assuming the Salpctcr initial mass function with the 
tabulated mass range. We also assume the escape fraction /esc = 
0.1 and the clumping factor C = 1 for the star formation density. 



Table 2 

Characteristic redshifts of NIR filters 



filter 


2^Lya;,in 


2Lya,out 


/ 


5.0 


6.6 


J 


8.0 


10.8 


H 


11.4 


14.0 


K 


15.4 


19.4 


L 


24.6 


30.6 



Note. ^Lyo.in and 2:Lya,out 

arc the redshifts at which the 
Lya break enters and goes out 
of the filter transmission, respec- 
tively. 



Table 3 

NIR afterglow intrinsic colors 



color 




(X v-'"/^ (p = 2.5) 


I - J 


0.70 


1.0 


I - H 


1.3 


1.9 


I - K 


2.0 


2.8 


I - L 


3.1 


4.3 


J - H 


0.62 


0.86 


J - K 


1.3 


2.8 


J - L 


2.4 


3.3 


H - K 


0.67 


0.91 


H - L 


1.8 


2.4 


K - L 


1.1 


1.5 



Note. — In the observer's frame, the NIR 
afterglow spectrum of the GRBs at 2: ~ 10 
is proportional to between ^ 1 minute 

and several hours after the initial burst and 
proportional to for later time. 



